Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



The central region of M31 observed with XMM-Newton* 
II. Variability of the individual sources. 

J. P. Osborne 1 , K.N. Borozdin 2 , S.P. Trudolyubov 2 , W.C. Priedhorsky 2 , R. Soria 3 , R. Shirey 4 , C. Hayter 1 , 
N. La Palombara 6 , K. Mason 3 , S. Molendi 6 , F. Paerels 7 , W. Pietsch 8 , A.M. Read 8 , A. Tiengo 5 , 

M.G. Watson 1 , and R.G. West 1 

1 Department of Physics & Astronomy, University of Leicester, Leicester LEI 7RH, UK 

2 NIS Division, Los Alamos National Laboratory, Los Alamos, NM 87545, USA 

3 Mullard Space Science Laboratory, University College London, Holmbury St. Mary, Dorking, RH5 6NT, UK 

4 Department of Physics, University of California, Santa Barbara, Santa Barbara, CA 93106, USA 

5 XMM-Newton SOC, VILSPA-ESA, Apartado 50727, 28080 Madrid, Spain 

6 Istituto di Fisica Cosmica "G.Occhialini" , Via Bassini 15, 20133, Milano, Italy 

7 Columbia Astrophysics Laboratory, Columbia University, New York, NY 10027, USA 

8 Max Planck Institut fur Extraterrestrische Physik, Giessenbachstrafie, D-85741 Garching bei Miinchen, 
Germany 

Received [date]; accepted [date] 

Abstract. We present the results of a study of the variability of X-ray sources in the central 30' of the nearby 
Andromeda Galaxy (M31) based on XMM-Newton Performance Verification observations. Two observations of 
this field, with a total exposure time of about 50 ks, were performed in June and December of 2000. We found 116 
sources brighter than a limiting luminosity of 6xl0 35 erg s _1 (0.3-12 keV, d = 760 kpc). For the ~60 brightest 
sources, we searched for periodic and non-periodic variability; at least 15% of these sources appear to be variable 
on a time scale of several months. We discovered a new bright transient source ~2.9' from the nucleus in the 
June observation; this source faded significantly and was no longer detected in December. The behaviour of the 
object is similar to a handful of Galactic LMXB transients, most of which are supposed to harbor black holes. 
We detected pulsations with a period of ^865 s from a source with a supersoft spectrum. The flux of this source 
decreased significantly between the two XMM observations. The detected period is unusually short and points 
to a rapidly spinning magnetized white dwarf. The high luminosity and transient nature of the source suggest its 
possible identification with classical or symbiotic nova, some of which were observed earlier as supersoft sources. 
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1. Introduction 



At a distance of 760 kpc ( van den Bcrgh 2000 : adopted 
throughout this paper), M31 is close enough to allow 
detailed study of individual sources within the galaxy 
using modern X-ray telescopes. Over 100 discrete X- 
ray sources in M31 were detected with the Einstein 



observatory (Trinchieri & Fabbiano 1991, hereafter TF; 
van Speybroeck et al. 1979] ). Primini, Forman & Jones 



(1993, hereafter PFJ) reported on the detection of 86 X- 
ray sources in the central 34' of M31 with the ROSAT 



HRI. Supper et al. (2001, hereafter 3u01) recently pub- 
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lished the complete catalog of 560 sources detected in 
a ~10.7 deg 2 survey of M31 with the ROSAT /PSPCQ 
From the extrapolation of the luminosity distribution, 
PFJ concluded that the detected population of X-ray 
sources could account for only ~15-26% of the unre- 
solved X-ray emission in M31, suggesting that the re- 
maining emission is truly diffuse or due to a new class 
of X-ray sources. Spectral analysis of ROSAT data per- 
formed by Borozdin & Priedhorsky (2000) and observa- 
tions with X MM (Shirey et al. 2 001, hereafter [Paper I| ) 
and Chandra (Garcia et al. 2001) showed that the unre- 
solved X-ray emission in the bulge of M31 is significantly 
softer than most of the point sources and can be approxi- 



This is an extension of the p reviou sly published catalog of 
396 sources (Supper et al. 1997, |Su97|). 



2 



Osborne et al.: M31 observed with XMM-Newton 



II. 



Variability 



o 



e 



5 M 



4 



30 



O 



Q 



0:42:00 



Fig. 1. The 0.3 - 10 keV EPIC PN 
image of the central region of M31. 
The data are from the 31-ks expo- 
sure taken on June 25, 2000. The 
intensity of the color is proportional 
to the logarithm of the number of 
counts collected (with a maximum 
of 5000 counts). PN chip bound- 
aries are clearly visible. The posi- 
tions of ten variable sources are en- 
circled and labelled with the num- 
bers listed in Table @. 



Table 1. XMM-Newton PV observations of the core of 
M31. 



performing spectral and timing studies in a field of point 
sources and extended emission. In Paper I we focused on 



Date & Time (UTC) 


Rev 


Obs. ID 


Exp.(ks) 


25/6/2000 (10:44-20:25) 
28/12/2000 (0:10-3:34) 


100 
193 


0112570401 
0112570601 


34.8° /30.7 6 
12.2 a /9.8 6 



a for EPIC-MOS 
6 for EPIC-PN 



mated with an optically-thin thermal plasma model (kT ~ 
0.3 keV) as expected from truly diffuse emission. 

In the first Chandra observation of M31, the nuclear 
source seen with Einstein and ROSAT was resolved into 
five sources (Garcia et al. 2000, hereafter p2000 ). One of 
these sources is located within 1" of the radio nucleus of 
M31 and exhibits an unusually soft X-ray spectrum, sug- 
gesting that it may be associated with the central super- 
massive black hole. A few more pairs of previously un- 
resolved sources and a new transient were also detected 
within 30" of the nucleus. 

We report on observations of the Andromeda 
Galaxy (M31) carried out with XMM-Newton 



(Jansen et al. 2001) during its Performance Verification 
(PV) phase. M31 was selected as an XMM PV target in 
order to demonstrate the capabilities of the mission in 



the group properties of the X-ray point sources and on 
the diffuse emission. The spectral properties of discrete 
X-ray sources in the XMM-Newton exposures will be 
discussed in Trudolyubov et al. (2001, Paper III). In 



this present letter we discuss the variability of individual 
X-ray sources in M31. In section |2|, we summarize the 
XMM-Newton observations and our data-reduction pro- 
cess. In section [| we discuss our detections of transient 
and periodic variability in individual sources. Finally, we 
present our conclusions. 

2. XMM-Newton Observations 

The bulge of M31 was observed with XMM-Newton on 
June 25, 2000 and again on Dec 28, 2000 (see Table Q). 
The observations were centered on the core of M31 (a — 
00 h 42 m 43!0, 6 = +41°15'46.0" J2000), with a field of 
view of 30' in diameter for the three European Photon 
Imaging Camera (E PIC) instruments . The two EPIC 
MPS instruments ([Turner et al. 200l|) and the EPIC 
PN ( Striider et al. 2001 ) operated in full- window mode 
with the medium optical blocking filter. The Optical/UV 
Monitor Telescope (OM; |Mason et al. 200l| ) filter wheel 
was set to the blocked position during the June observa- 
tion. During the December observation, exposures were 
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# 


Source name a 


a (2000) 


5(2000) 


.Tun 


r Dec 


Comments 11 


1 

1 




nn-/t9-n c i q 


41 


13:29 


<^9 7 


1 9 ft-|-9 A 
1Z.UHZZ.U 




o 
A 




nn-/i9-i 9 ^ 


41 


18:00 


l.oznl.o 


1 7 ^4-/1 1 




Q 

o 


XA/TA/TTT T0049' : 14 1 4-41 1 SDS 


UU.4:i.04:. J. 


41 


18:08 


07 q+i 7 


^fi 4 




4 


CXO 1004242 0+41 1 60S 


00-42-42 2 


41 


16:09 


58.1±2.1 


45 fi+3 


Chandra transient ("02000 and S3 2i 


5 


RX J0042. 7+4116 


00:42:47.2 


41 


16:28 


9.5±1.3 


94.0±3.6 


TF#59, PFJ#50, Su97#198, 
Su01#195 


6 


XMMU J004247.5+411158 


00:42:47.5 


41 


11:58 


4.0±1.3 


<3.0 




7 


CXOU J004257. 1+41 1843 


00:42:57.1 


41 


18:43 


<2.0 


6.3±2.2 


seen on Oct 13, 1999 with Chandra 


8 


RX J0043.3+4120 


00:43:18.9 


41 


20:19 


23.4±1.3 e 


2.5±1.4 e 


SSS / , TF#87, Su97#235, Su01#235 


9 


XMMU J004319.4+411759 


00:43:19.4 


41 


17:59 


43.6±1.7 e ' 9 


<5.r- h 


SSS 7 , 865-s pulsations (§3.3) 


10 


RX J0043.4+4118 


00:43:27.9 


41 


18:35 


8.6±0.8 e 


6.5±1.5 e 


SSS 7 , TF#89, PFJ#80, Su97#240, 



Su01#241, SNR 



a The names for ROSAT and Einstein sources were adopted from SuOl. If the source had not been found in SuOl, we 
used the acronym XR followed with truncated XMM coordinates. The sources seen first with Chandra or XMM are named 
according to the naming conventions for these missions. 
6 Source count rates for the XMM-Newton observations in June and December of 2000 are given in units of 10 _3 counts/s 
(MOS1) and are corrected for vignetting. Quoted count rates are measured in the 0.3-10 keV energy range, except for 
supersoft sources (SSS), which are measured in the 0.3-1.5 keV band. In cases when sources are not detected, 2a upper 
limits are listed. 
c References to the earlier observations: 



TFj = Trinchieri fc Fab biano (1991); |PFJ 
^u97| = Supper et al. (1997); |Su0l] = Supper et al. (2001); |G2000| = Garcia et al 



= Primini, Forman & Jones (1993); 



(2000). 

Source near the chip edge; the detected count rate might be an underestimation of real count rate. 
Source count rates for 0.3-1.5 keV band. 
Supersoft source. 

With contribution from a nearby faint source. 

Upper limit is defined by the contribution of a nearby source. 



obtained with the B and UVW1 filter of the OM; results 
of these observations will be presented elsewhere. 

A background flare occurred during the final 5 ks of 
the EPIC exposures from the June observation. Data ob- 
tained during this background flare were excluded from 
the variability studies. 

We used the XMM-Newton Science Analysis System 
(versions 4.1 and 5.0.1) to reduce the EPIC data to cal- 
ibrated event lists, produce images, and extract light 
curves. A combination of SAS programs and external soft- 
ware was applied to further analyze the data. The data 
from MOS1, MOS2 and PN were extensively compared to 
confirm the consistency of our results. The variability of 
individual sources within the longest (June) observation 
was studied by Fourier analysis. The value of the pulsa- 
tion period detected in one of the supersoft sources (see 
below) was refined by the epoch-folding method. 



3. Individual X-ray sources 

The EPIC images of the central 30' of M31 (Fig. | rep- 
resents part of EPIC PN image) contain more than 100 
discrete X-ray sources as well as unresolved emission near 
the centre (see Paper I for more details). The detections 
include sources seen with Einstein, ROSAT, or Chandra 
( fff[ |PFJ| , |Su97j , j5u0l| , |G2000| ) as well as ne w sourc es. 
The transient source discovered by Chandra ( |G2000 ) is 



clearly seen in both XMM-Newton observations (§3.2). 
In the June observation, there is a bright new transient 
which was not seen previously (§3.1). Several other bright 
sources demonstrated significant variability on half of a 
year time scale, so that uncertainties of count rate mea- 
surements do not overlap for the two observations. We 
have analyzed 66 sources with limiting count rate^ 4 
counts/ks and found 10 variables. These sources are listed 
in Table || along with their XMM positions, June and 
December count rates in MOS1, and any relevant iden- 
tifications (though only count rates in MOS1 are cited, 
we selected the sources based on the data from all three 
XMM-Newton instruments). The source positions were 
derived from comparison with Chandra calibrated images 
of the central region of M31. The estimated positional er- 
ror of 1"— 3" (in diameter) depends on the distance of the 
source from the boresight axis, its intensity and nearby 
sources. We estimate 3" as a conservative error limit for 
the XMM positions presented in this paper. 

Taking into account the total number of sufficiently 
bright detected sources, we conclude that at least ~15% 
of all sources in the central field of M31 are variable on 
a time scale of several months. We believe, however, that 
this value should be considered a lower limit, because the 



2 Some of the sufficiently bright sources were excluded from 
our analysis because of the confusion with nearby sources or 
because they were too close to the chip edge 
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sensitivity of our analysis was a function of the source flux, 
and the variability of many of the fainter sources would 
not be detected. Below we discuss in more detail several 
individual variable sources in M31. 



3.1. X-ray Nova in M31 

We found a new bright X-ray source in EPIC PN and 
EPIC MOS images from the June observation. The co- 
ordinates of the new source are presented in Table ^|. 
The source was neither detected in previous Einstein 
and ROSAT observations nor has been reported from 
Chandra observations. During the June observations the 
flux from the source was (1.57±0.09) x 1CP 13 erg s _1 cm~ 2 
(0.3-10 keV), which corresponds to a luminosity of ~ 
1.1 xlO 37 erg s _1 in the c ited energ y band assuming 
a power-law spectrum (see Paper III ). The source flux 
faded to below the background level before the next 
XMM-Newton observation of the same field on Dec 28, 
2000. The upper limit for the count rate for this second 
observation was less than 6.4 counts/ks (2er upper limit for 
MOS1). The detected luminosity of the source during the 
June 25 observation, as well as its several-month time scale 
to fade to quiescence, is typical for bright X-ray transients 
in our Galaxy (see reviews by Tanaka & Shibazaki 1996 



and Chen, Shrader, & Livio 1997) 



3.2. Chandra transient 

A bright transient source designated 
CXO J004242.0+411608 was discovered with Chandra 
during a 17.5 ks observation of the core of M31 
on Oct 13, 1999 ( |G2000| ). We detected an X- 
ray source at the same location during both our 
XMM-Newton PV observations. The flux of the 
source measured on June 25, 2000 was equal to 
(5.5±0.5)xl0~ 13 erg s" 1 cm~ 2 (L x ~ 3.8xl0 37 erg s" 1 ). 
The source was detected again on Dec 28, 2000 with flux 
of (4.2±1.0)xl0 -13 erg s _1 cm~ 2 . Including both the 
Chandra and XMM detections, the source has remained 
bright for more than 14 months, with only slight possible 
fading (^25%) between the two XMM-Newton observa- 
tions separated by half of a year. The X-ray luminosity 
and spectral shape detected by XMM-Newton ( Paper III ) 
correspond to those of the har d/low state of Galacti c 
black hole transients (e.g., Tanaka fc Lcwin 1995|) . 
Typically Galactic X-ray tr ansients fade significantly in 
less than ~100 days (see Chen, Shrader, fc Livio 1997 
for a review); however, a long plateau in a hard spectral 
state was observed from t h e X-ray transient GRS 1716- 
249 dSunyaev et al. 1994 |Revnivtsev et al. 1998| ). We 
should mention that our sampling does not allow us 
to distinguish reliably between single and repeated 
outbursts. It will be interesting to follow the evolu- 
tion of CXO J004242. 0+411608 during future planned 
observations with XMM-Newton . 




0.0001 0.0010 0.0100 

Frequency (Hz) 

Fig. 2. Power density spectra for a supersoft source in 
M31 obtained from each of the three EPIC instruments 
during the observation on June 25, 2000. The peak power 
from each of the three instruments occurs at a frequency 
corresponding to to a period of ^865 s. The vertical dotted 
line denotes the best fit value of the period found by the 
epoch-folding method. Detection thresholds are indicated 
by horizontal dashed lines labelled with the approximate 
probability that any of the 1500 or 1300 frequency bins in 
each PDS (for MOSs and PN respectively; no averaging 
applied) wou ld have a noise value exceeding the indicated 
power level ( Vaughan et al. 1994). 



3.3. Supersoft 865-s pulsator 

We detected significant oscillations in the X-ray flux from 
one of the bright variable sources in our field (source #9 in 
Table |). The Fourier power density spectra (PDS) from 
each of the three independent EPIC instruments (Fig. |2|) 
show a highly significant periodic signal from this source 
during the June observation, with a best-fit period of 
865.5±0.5 s. The folded light curve for this source (Fig. |||) 
is quasi-sinusoidal with an amplitude of about 40%. 

In June 2000 this 865-s pulsator had a spectrum similar 
to the three sources in the XMM field previously identified 
by Kahabka (1999) as supersoft source (SSS) candidates 
based on their ROSAT /PSPC hardness ratios and lumi- 
nosity (one was identified with a supernova remnant). The 
temperature of the blackbody spectral fit was kT(,f,=61±2 
eV with reduced x 2 =1.5 (for EPIC-MOS data). All four 
SSS exhibit blackbody-like spectra with effective tempera- 
tures in the range kT~50-150 eV and no detectable X-ray 
emission above ~1.5 keV (see |Paper III for more detailed 
discussion of spectral parameters). Such sources are typi- 
cally interpreted as accreting white dwarfs in binary sys- 
tems, powered by nuclear burning of the accreted matter 
on their surfaces (see Kahabka fc van den Hcuvcl 1997 for 
a review) . 

The source count rate during the June 25 observation 
was equal to 44 counts/ks in the 0.3-1.5 keV energy range 
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(MOS1, see Table|^). By Dec. 28 the source count rate had 
faded down to a level below 5 counts/ks (MOS1); thus, it 
was not possible to detect pulsations during the second ob- 
servation. All three supersoft sources in the NE portion of 
our field (sources 8, 9 and 10) faded between the June 
and December observations (see Table ||). We are aware of 
a possible sensitivity degradation in this region of MOS1; 
however, the result is confirmed by MOS2 and PN data. In 
particular, the count rate of the supersoft pulsator as de- 
tected by MOS2 dropped down from 41±1.5 counts/ks in 
June to 3±1.5 counts/ks in December, while many other 
sources in the field showed no significant variability or be- 
came brighter in December. 

We note that there is a nearby source (a = 
00 h 43 m 21!2, 6 = +41°17'52", J2000) with a harder spec- 
trum which is much fainter than the supersoft pulsator in 
the June observation but dominates this region of the sky 
for the December observation. The distance between two 
sources is about 10", so that they are spatially resolved 
with XMM-Newton , but their counts are not completely 
spatially separated on the detector. Because the relative 
flux of the harder source was much lower during the June 
observation, it did not significantly affect our timing and 
spectral analysis of the SSS. The ROSAT /PSPC would 
barely be able to resolv e the two sources, so RO SAT's 
RX J0043. 3+4117 flSu0l| #236) may include contributions 
from both of them, although its position coincides with 
that of the harder source rather than the SSS pulsator. 

The 865. 5-s period is the shortest among all known 
SSSs (Greiner 2000). Interpreted as the binary orbital 
period, it would be too short to accommodate a main- 
sequence companion and would suggest a degenerate sec- 
ondary. It may be more plausible to assume that the pul- 
sations indicate that the white dwarf possesses a mag- 
netic field large enough to modulate the X-ray emission 
yet not so large that the spin and orbital periods are 
locked, e.g., as in intermedia te polars or DQ Her stars 
(see review by Cordova 1995| ). However, the luminosity 
of the object ~ 1.7xl0 37 erg s^ 1 (0.3-1.5 keV) is sev- 
eral orders of magnitude higher then typical for inter- 
mediate polars luminosity range Lx = 10 31 -10 (see e.g. 
Patterson 1994). The high luminosity and transient na- 
ture of the pulsator may indicate steady burning in a 
post-nova stage, as has been observed in a few classical 
and symbiotic novae (Kahabka & van den Hcuvcl 1997). 



An alternative explanation for the nature of the SSS could 
be a double degenerate polar similar to RX J1914. 4+2456 
( |Haberl fc Motch 1994 [Ramsay et al. 200C| ). An interpre- 
tation of the source as a foreground object is unlikely due 
to the lack of an optical counterpart in XMM OM images 
up to the limiting magnitude of ~19 in B filter. 



3.4. Variability of other sources within a single 
observation 

We have searched for coherent periodic modulation on 
time scales from ^10 to ~1000 s for about 60 of the bright- 




Phase (TJD 1 1720 = 0) 

Fig. 3. Light curve of the pulsating supersoft source in 
M31, folded using a period of 865.5 s. The data is from 
the June 25, 2000 observation, and counts from all three 
EPIC instruments have been used. Contribution of the 
background is within the error bars. 



est sources in the XMM-Newton field of view; however, 
only in one case (see previous subsection) was significant 
variability detected. The 90% upper limits to periodic 
modulation fractions, calculated as the ratio of the sine 
amplitude to the constant flux level for periods of 10 000, 
300, and 10 s were obtained in all other cases. These upper 
limits vary from 3.6% for the brightest source to ~30% for 
the faintest sources of the sample. 

The lack of detectable variability for many of the in- 
dividual sources in M31 may seem surprising compared 
with rich variability observed from Galactic sources. Our 
observations so far have included mainly the bulge of M31, 
where low-mass X-ray binaries (LMXBs) are most preva- 
lent. Such systems commonly show dips, bursts, and quasi- 
periodic variability rather than coherent pulsations. We 
expect to detect more X-ray pulsars during planned XMM 
observations of M31 in fields along the disk of M31, where 
population I stars and high-mass X-ray binaries (HMXBs) 
dominate. We must note also that the sensitivity to vari- 
ability depends strongly on the brightness of the source, 
and hence our data are not very sensitive to the variabil- 
ity of the many faint sources. Finally, we did not perform 
timing analysis for time scales shorter than ~10 s because 
the absolute time calibration was not possible with the 
Observation Data Files available. 

We have also looked for non-periodic variations on all 
accessible timescales for the same set of objects. No X-ray 
burst has been detected, however our sensitivity to bursts 
is restricted to a relatively small luminosity interval by 
the low count rate of most sources and by the Eddington 
limit. 
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4. Conclusions 

In this letter we present results obtained from 
XMM-Newton PV observations of M31. Two observations 
separated by half of a year were carried out in June and 
December of 2000. This letter is the second paper in the 
series describing results of these observations. 

Significant variability of individual sources was de- 
tected both between the two XMM observations and in 
comparison with earlier results of other missions. At least 
~15% of the sources appear to be variable, and we con- 
sider this value to be a conservative lower limit. 

A new bright transient source was detected during the 
observation of June 25 but faded before the December ob- 
servation. Probably it is an LMXB transient source, sim- 
ilar to a handful of such sources observed in our Galaxy, 
most of which are supposed to harbour black holes. 

Another transient source, first detected by Chandra 
( G2000 ), was bright during both XMM-Newton observa- 
tions. The flux of the source did not change significantly 
in observations separated by six months, which is not 
typical for Galactic X-ray transients but is reminiscent 
of the behavior of the black hole candidate GRS 1716- 



249 during its 1993-19 94 outbursts flSunyaev et al. 1994 , 
Revnivtsev et al. 1998). 

X-ray pulsations with period of ~865.5 s and quasi- 
sinusoidal pulse profile were detected from one of the su- 
persoft sources in our field. It was bright during the June 
observation, but had faded such that it became unde- 
tectable in December. The period of the detected pulsa- 
tions is the shortest among known SSSs. A likely source 
of the pulsations is a magnetized rapidly spinning white 
dwarf; however, the luminosity of the source is much 
higher than for typical CV systems. The detected X-ray 
flux may be generated by steady nuclear burning in a post- 
nova stage of a classical or symbiotic nova. 
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